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Abstract. Using a classical ensemble method, nonsequential double ionization is 
predicted to exist with elliptical and circular polarization. RecoUision is found to be 
the underlying mechanism and it is only possible via elliptical trajectories. 
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1. Introduction 

An anomalously high degree of electron correlation is well known to be present in 
multiphoton atomic double ionization (see reviews [T] [2]). Two ionization channels 
are available: the atom may either lose the two electrons one by one, which is called 
sequential double ionization (SDI), or lose the two electrons at the same time, which 
is called nonsequential double ionization (NSDI). Compared to SDI, which can be 
well explained by the ADK tunneling theory j3] based on a single active electron 
approximation, NSDI is a much more complicated (and thus interesting) phenomenon 
of atomic physics. Recently, open questions in this domain associated with departures 
from the usual linear polarization of the laser pulses have received close theoretical 
attention [4l,l5l,l6]. 

The accepted scenario for such high 2e correlation is a three-step model in which 
one electron is ionized via tunneling when the Coulomb barrier is tipped down by the 
laser electric field. Second, the electron is carried away and then driven back to the ion 
core as the laser field reverses phase. Third, the energized returning electron recombines 
at the nucleus to generate high harmonics ^ , or kicks out a second electron [8j to 
provide double ionization, as sketched in Fig. [TJ 

This classically-inspired picture leads to the understanding that little or no NSDI 
can be expected if the laser field is elliptically polarized, since the first electron 
cannot return to the ion core under elliptical polarization. This is at odds with some 
experiments. Characteristic NSDI events have been observed with the molecules NO 
and O2 [9j and atomic magnesium [10] under circular polarization. We report here 
an explanation for these phenomena, and show why they do not require invention of 
a diflFerent non-collision scenario. To be clear that we are dealing with NSDI as it is 
usually understood, we show in Fig. [2] that the characteristic knee signature for NSDI 
is present for all degrees of ellipticity, up to and including circular polarization. 




Figure 1. High harmonic generation and double ionization via the three-step scenario, 
in which the reversed laser phase directs the first ionized electron back to the core. 
Linear polarization of the laser pulse is imphed, in order that the returning electron 
collide with the nucleus (see \1^S^). 
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Figure 2. Ionization probability obtained in numerical experiments as a function of 
laser peak intensity under different ellipticities for a model "Xe/Kr" atom. So-called 
knee structure is evident for a full range of ellipticities. 



2. Elliptical trajectories 

To understand NSDI with elliptical polarization we have undertaken a large series of 
numerical "experiments" (see last section below for details), in which two electrons 
are followed along classical trajectories. Successful NSDI events were back-analyzed to 
determine the sequence of events that lead to an NSDI conclusion. Back-analysis of 
classical trajectories is a very powerful tool to investigate detailed physical processes 
fn\ and a quantum analog using wave function masking has also been developed by 
Haan, et al. [12]. 
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Figure 3. Two typical recollision trajectories for ellipticity 0.5. The blue boxes 
of the left two figures are greatly expanded to give the right two figures. Top: the 
recolliding electron's trajectory is an almost perfect ellipse. It returns to the origin in 
one laser cycle and kicks out the second electron. Bottom: the elliptical trajectory of 
the recolliding electron is distorted by the ion's Coulomb force and executes several 
cycles before coming close enough to the ion core to kick out the second electron. The 
arrows and numbers are used to show the temporal motion of the first electron. 
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The surprising conclusion of this back analysis is that every successful NSDI 
trajectory is elliptical! That is, in all cases examined, the first ionized electron circulates 
along an elliptical trajectory for one or more cycles and then collides and ejects the 
second electron. Typical trajectories showing one-cycle and multi-cycle ellipses are 
given in Fig. [3) Thus the answer to the question raised by the apparent failure of 
the recollision scenario for polarizations that are not linear is that there is no failure 
- recollision is still the mechanism. More than that, recollision is possible only via 
elliptical trajectories. 



3. Understanding the elliptical trajectories 

The dominance of the elliptical trajectories can be understood from a simple analysis 
based on free electron kinematics. Our analysis begins with the common assumption 
that electron motion after ionization is adequately described by free electron physics in 
the presence of the laser field. The core's Coulomb field plays a role, but it is secondary, 
so the force felt by an electron comes mainly from the elliptically polarized field 

E{t) = Eof{t)[x8m{ujt + (p)+£yco8{ujt + (/))], (1) 

as sketched in Fig. [4} Here is a random phase from shot to shot in our series of 
numerical laser- atom interactions, and f{t) is a slowly varying several-cycle envelope 
whose time dependence can be ignored in first approximation. 




Figure 4. Illustration of an elliptically polarized laser pulse. The red curve is the 
electric field along the x direction and the blue curve is the electric field along the y 
direction. The ellipticity is 0.5 here. 



The charge of the electron is -1 in atomic units so the free-electron velocities along 
the X and y directions are (having put / = 1): 

rt 



and 



Ho 

Vox H — ^o[cos(cL;t + (/))- cos(cL;to + 0)], (2) 

id 

Vy{t) = voy EQ[sm{(jjt + 0) - sm{ujto + 0)]. (3) 

The first terms are the initial velocities at the time of ionization, independent of the laser 
force. By the nature of first escape over the tipped Coulomb barrier (or via tunneling 
through it), the initial electron velocity in the polarization direction is zero. So Vqx = 0. 
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However, the initial velocity in the transverse direction can be nonzero. The exact 
form of the transverse distribution is not important, and has only been derived in a 
one-electron model [13], but its existence is important and has begun to be taken into 
account (see a recent treatment of ionization under elliptical polarization by Shvetsov- 
Shilovski, et al. [4J). The second terms are the electron oscillation velocities in the field 
and the last terms are laser-induced drift velocities. 

Because the ionization is strongly dependent on field strength the first electron is 
most probably ionized at the peak of the field strength. To a first approximation, for 
ellipticities not close to circular, this makes sin(cL;to + 0) ^ ±1 and cos(cL;to + 0) ^ 0. 
Under this scenario the electron has the drift velocity Vyd in the y direction: 

Vyd = VOy ± . (4) 

(jU 

This drift velocity has been used to obtain the first formula for NSDI probability as a 
function of ellipticity [14j. Apart from this drift motion, one sees that the trajectory is 
a perfect ellipse, which we call the "circulation ellipse", with half major axis Eq/uj'^ and 
half minor axis sEq/uj'^. 

For intensity 6 x lO-*^^ W/cm^ and wave length 780 nm, we have Eq/uj'^ ^ 40 a.u. 
so the field-dependent part of the transverse drift distance for one laser cycle is about 
(sEq/u)) X [2ti/uj) = 27r£ Eq/u)'^^ which is about 6 times larger than the size of the 
circulation ellipse. This means that the electron will quickly be taken away from the 
ion core and no recollision could occur. However, from the range of velocities Voy that 
are available, there will be a compensating value, allowing a perfect elliptical oscillation 
that returns the first electron to the ion core, for a collision that may kick out the second 
electron. 

4. The role of the ion core Coulomb potential 

The ion core Coulomb potential, which is neglected in the free electron analysis of the 
previous section, plays a role in increasing the probability of NSDI. As explained above, 
balancing the laser-induced drift velocity along the y direction with a compensating 
transverse initial ionization velocity is the reason for successful recollisions. Exact 
balancing yields a purely elliptical trajectory, and a nearly perfect example is shown 
in Fig. |3] (a) and (b). If the balancing is not so precise, the first ionized electron will 
drift away, although the drifting velocity is small. With the aid of the ion core Coulomb 
potential, however, slow drifts will be pulled back. Slowly drifting circulation ellipses 
can be pulled back and forth around the ion core several times before the first electron 
finds the second one and kicks it out. This is shown in Fig. |3] (c) and (d). In fact, we 
find in our numerical experiments that most of the successful recollision trajectories fall 
into this category. 
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5. Method 



The method that we have used to obtain the numerical experiment data presented is the 
classical ensemble method, which has been described in detail in [15j. The contribution 
of the classical ensemble method to interpretation of double ionization phenomena is 
well known and has been critically analyzed p!6l [TTj . A microcanonical but otherwise 
essentially random phase-space-filling distribution of initial conditions (typically 10^-10'' 
in number) is generated using a many-pilot-atom method [18j before turning on the laser 
field. 

The classical atom model has a Z = +2 charged ion core fixed at the origin, and 
the energy of the two-electron bound system is given by 

^tot = y + y - ^Kc(ri, a) - ZVUr2, a) + VUn - rl, 6), (5) 
with Vsc standing for the soft-core Coulomb potential of the familiar "Rochester" form 

ysc{r,c) = ^^== (6) 

The energy of each 2e member of the ensemble is set to be -1.3 a.u., which is close to 
the binding energies of both Kr (-1.41 a.u.) and Xe (-1.23 a.u.). (The binding energy of 
a 2e system is calculated as the negative sum of the 1st and 2nd ionization potentials. 
For example, the 1st and 2nd ionization potentials of Kr are 13.999 eV and 24.359 eV, 
which yields a binding energy of -(13.999+24.359) eV = -38.358 eV = -1.41 a.u.) 

An elliptically polarized pulse as shown in Fig. |4] is then applied to the system 
and at each moment the positions and momenta of both electrons are recorded, as 
determined by solution of the relevant coupled Newtonian equations. Experience with 
this method [20] has shown that out-of-plane eflFects can be neglected under current 
experimental conditions so we only work on the x-y plane. 



6. Summary 

In summary, we have shown that the characteristic NSDI knee structure still appears 
with elliptical and circular polarization. This is consistent with results of highly 
cited experiments 00^. Elliptical polarization leads to elliptical trajectories and our 
numerical experiments show that all successful NSDI events come from recollisions at 
the end of an elliptical trajectory. The dynamical process determining these recoUiding 
elliptical trajectories is explained using free electron kinematics. The residual effect of 
Coulomb attraction to the nucleus plays a secondary role in many cases. 
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